Distribution of spectral weight in a system with disordered stripes 
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The "band-structure" of a disordered stripe array is computed and compared, at a qualitative 
level, to angle resolved photoemission experiments on the cuprate high temperature superconductors. 
The low-energy states are found to be strongly localized transverse to the stripe direction, so the 
electron dynamics is strictly one-dimensional (along the stripe). Despite this, aspects of the two 
dimensional band-structure Fermi surface are still vividly apparent. 



I. INTRODUCTION 



There is strong evidence, in part based on careful anal- 
ysis of single particle spectral functions deduced from 
angle resolved photoemission (ARPES), that the nor- 
mal state of the high temperature superconductors is not 
a conventional Fermi liquid. At the same time, there 
is clear evidence that in some, and possibly in all of 
these materials, there are local, self-organized quasi-one- 
dimensional structures, "stripes" , which substantially af- 
fect the electron dynamics. Many striking features of the 
ARPES spectrum have a natural interpretation in terms 
of an electronic structure dominated by stripes [p]-p|. 

One such feature, which is reproduced in Fig. El, is that 
the low-energy spectral weight lies predominantly in re- 
gions of the Brillouin zone (BZ) in the vicinity of the 
(±7r, 0) and (0, ±7r) points. (These are often referred to 
as the "anti-nodal regions" as this is where the maximum 
in the d-wave superconducting gap occurs.) These re- 
gions have very straight boundaries which run parallel to 
the momentum axes and are displaced by approximately 
±7r/4 from them. This spectral weight has been inter- 
preted in terms of quarter filled stripes (so hp = 7r/4) 
along the a and b directions of the lattice • 

Bowever, a number of features of the ARPES 
data are inconsistent with the most naive picture 
of an array of strictly one-dimensional quarter filled 
stripes. Of these, the most obvious is the occur- 
rence of features, especially in the "nodal regions" 
of the BZ [around (±7r/2,±7r/2)], which are reminis- 
cent of the underlying, two dimensional band struc- 
ture. An example is the appearance of low-energy spec- 
tral weight along segments parallel to the lines con- 
necting the anti-nodal points in I^-^Sr^CuO^LSCO) 
and Lai.6-xNdo.4Sra;Cu04(LNSCO). A similar signature 
also appears in Pb-doped Bi2Sr 2 CaCu208+<5(BSCCO), 
as shown in Fig. [I]. 

It is the principal purpose of this paper to elucidate 
how this apparently two-dimensional structure arises nat- 
urally from a slightly more sophisticated analysis of the 
band-structure of a stripy system. In doing so we ex- 
tend the initial work of Ref. ^ and complement a re- 
cent study of the effects of realistic band parameters on 
the spectrum by Fleck, Pavarini and Andersen M. We 



show that a system can exhibit what seems to be a Fermi 
surface of a two-dimensional metal despite the fact that 
the dynamics of its low-energy electrons is entirely one- 
dimensional ! This conclusion holds even in the presence 
of interactions that preclude any quasiparticle-like de- 
scription of the system. Under such conditions although 
the low-energy features in momentum space are sharp, 
the spectral function, considered as function of energy 
for fixed momentum can be broad km] 



Our results are easily summarized. We have computed 
the band-structure of electrons in the potential gener- 
ated by a typical configuration of the "slow" collective 
fields that define the stripe order. (The explicit Bamil- 
tonian is presented in Section |l[) Figs. || and || depict 
the k space distribution of low-energy spectral intensity 
in the first BZ for disordered and ordered arrays of par- 
allel stripes, respectively. In both cases most of the low- 
energy spectral weight is concentrated in the anti-nodal 
regions. Bowever, while the ordered stripe array exhibits 
a spectral gap in the nodal region, the disordered array 
has low-energy weight there, much like that seen in ex- 
periment. 

We picture the disordered array as being a "snapshot" 
of a "fluctuating stripe array," although it could also 
reflect the effect of quenched disorder. The low-energy 
electronic states that are responsible for its spectral map, 
shown in Fig. ||, lie within the Mott gap. They are bound 
states that decay exponentially in the direction perpen- 
dicular to the stripes. Bowever, they are extended along 
the stripes, so the low energy electron dynamics is strictly 
one-dimensional. The resulting band-structure, which is 
shown in Fig. g, is, to a good approximation, a superpo- 
sition of the band-structures of single isolated stripes and 
of small clusters (typically pairs) of proximate stripes. 

The band-structure and spectral distribution of a sin- 
gle stripe are shown in Figs. || and [7| Bere, the fact 
that the electron wavefunctions in the gap have a non- 
negligible extent transverse to the stripe can be seen to 
produce an image of the full two-dimensional band struc- 
ture. Bowever, as is apparent from the spectrum, unless 
an isolated stripe is nearly half-filled (which is physically 
implausible), it will not have any low-energy spectral 
weight near (tt/2, tt/2). This is the reason for the absence 
of spectral weight in the nodal region of the ordered ar- 
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ray since its wavefunctions are Bloch states constructed 
from single stripe bound states. 

In the disordered array of stripes, the spectral weight in 
the nodal regions originates primarily from anomalously 
close pairs of stripes. This is demonstrated in Figs. ||- 
[Toj , which show that for an isolated bi-stripe, there is 
low-energy spectral weight in the vicinity of (7r/2, it/2), 
even when the stripes are roughly quarter filled. The 
idea that the intensity of the nodal spectral weight is re- 
lated to the degree of stripe disorder jjj conforms with 
the experimental finding (see Fig. [I]) of more pronounced 
Fermi segments in the nodal direction in optimally doped 
LSCO, where the stripes do not statically order, than in 
LNSCO, where long-range stripe order is seen in neutron 
diffraction jlO|]. Several other experimental observations 
can be understood in terms of our results, as discussed 



in Sec. VI 



When residual interactions between electrons on a 
stripe are considered, each non-interacting localized band 
gives rise to a one-dimensional Luttingcr liquid. As in the 
non-interacting case the two-dimensional k space struc- 
ture, and especially the Fermi surface, remain prominent 
features of the low-energy spectral response. However, 
along a cut in the BZ perpendicular to the Fermi sur- 
face, the spectral function mimics Luttinger liquid be- 
havior. (See, e.g., the behavior along the nodal direc- 
tion in Fig. [□]. ) This justifies the application of one- 
dimensional physics along directions that are not neces- 
sarily aligned with the stripes, as was done recently in 
Ref. |. 

It is important to stress that we view the present re- 
sults as reliable for exploring the qualitative effects of 
stripes on electronic structure, but not as a realistic study 
of the cuprate superconductors. The parameters in our 
model have not been carefully adjusted to optimize any 
sort of fit to the data. We have certainly not included any 
"realistic" band-structure effects, such as second neigh- 
bor hopping, nor, except in Sec. |v|, considered any 
strong correlation effects, other than the stripes them- 
selves. Phonons, dynamical stripe fluctuations, effects of 
transverse deformations of the stripe potential, dynami- 
cal magnetic fluctuations, and all other forms of of static 
or dynamical disorder are neglected in our calculations. 
They will all certainly have important consequences for 
the details of the measured electronic structures. How- 
ever, it should also be stressed that, throughout this pa- 
per, we will be concerned with intermediate energies and 
length scales. Thus, for our purposes, the distinction 
between long-range and mesoscale stripe order, and be- 
tween static and fluctuating stripes is unimportant, al- 
though of course these distinctions are essential at low 
energies and long wave-lengths. 



II. THE MODEL 

Since the high temperature superconductors are doped 
antiferromagnets, it is reasonable to expect slow fluctua- 
tions of a collective field representing the local staggered 
magnetization. One can think of this field as resulting 
from a Hubbard-Stratonovich transformation of the in- 
teracting many-body problem. However, we do not solve 
a Hartree-Fock theory for this field, i.e., we do not find 
the configuration that minimizes the Hartree-Fock en- 
ergy. Instead we take as a minimal model a set of non- 
interacting electrons on a square lattice interacting with 
a static, staggered field which represents a characteristic 
"snapshot" of the field configuration: 

H = - X! ^ C x,yM C x+l,y,a + c x,y,a c x,y+l,a + H.C.) 
x,y,a 

+ ^(-l) ;E+y m(x, tf)4 jW ,«,Cr,y,«r , (1) 

x,y,a 

where c XjV ^ is the electron destruction operator at site 
(x, y) and spin a = ± and we have chosen units such 
that the hopping matrix element and the lattice constant 
equal 1. Specifically, we consider the simplest possible 
ansatz to represent stripe configurations in which 



m{x, y) = to J^J 0(x — x s ) 



(2) 



where to is a constant, {x s } are a given set of positions 
of anti-phase domain walls, and where O is the antisym- 
metric step-function: Q(x) — —Q(~x) = 1 for x > 
and O(0) = 0. This corresponds to an array of perfectly 
straight, site centered stripes of width one oriented in 
the y direction. Depending on the choice of the set {x s } 
the potential can either be regular or disordered in the 
transverse (x) direction. 

Some features of the solution below do depend on de- 
tails of this choice, such as whether the stripes are site 
or bond centered, whether they are of width 1 or wider, 
whether we include an additional collective field that cou- 
ples to the charge density 0], etc.; effects that we have 
explored, to some extent. However, the important qual- 
itative physics is apparent in this simplest of models, so 
we only report results for this model. Two more serious 
omissions, which are beyond the scope of the present pa- 
per, are the neglect of effects of the dynamical character 
of the collective fields, and the neglect of any shape de- 
formations of the stripe order. The latter approximation 
implies that the electronic states are Bloch waves in the 
y direction, with a wave- vector, k y , which is conserved 
mod 7r (due to the presence of the staggered field) . 

Without the domain walls the system has an energy 
gap of magnitude 2to and two bands with energies E = 
±\ / m 2 + e 2 , with e = — 2[cos(k x ) + cos(k y )]. An isolated 
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domain wall generates mid-gap states which are local- 
ized in the direction orthogonal to the wall, as is shown 
analytically, below. 

For an array of domain walls, it would be exceedingly 
difficult to find analytic solutions for the mid-gap states. 
Instead we diagonalize such systems numerically for a 
given realization of the model. In each case, we take 
m = 1 which corresponds to an "intermediate coupling" 
value of the ratio of the energy gap to the band-width 
of 2/8=1/4. We study systems with finite width, L x , in 
the x direction where typically L x = 320; according to 
all tests we have applied, this is large enough to elimi- 
nate most finite size effects. Since k y is a good quantum 
number, the only place the finite size in the y direction 
enters our calculations is when we perform sums over k y ; 
here, it is easy to show that the system sizes we have con- 
sidered, L y — 320, are more than adequate to eliminate 
finite size effects. 



III. BAND-STRUCTURE OF A STRIPE ARRAY 



In this section, we consider the spectral weight distri- 
bution and band structure of a stripe array. If one ig- 
nores details concerning the matrix elements describing 
the photo-excitation process one can infer directly the 
single-particle spectral function, A(k,u), from ARPES. 
For non-interacting electrons, A(k, ui) can be expressed 
in terms of the exact single particle energy eigenstates 
^/ a (r) and energies E a : 

A{k,u) = (L^-^e^-^i^rV) , 
A(f,f,cj) = Y f ^ a {i f )^ a (f)6(u-E a ) . (3) 

a 



In order to characterize the momentum distribution of 
the low-energy spectral weight and reveal features such 
as the Fermi surface, it is appropriate to integrate the 
spectral function over a narrow energy window below the 
chemical potential /i 

I(jfe) = [ A{k,u)duj . (4) 

J [i — Aio 

The experimental data in Fig. |l| and the theoretical re- 
sults in later figures are expressed in this way. 
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FIG. 1. Distribution of low-energy spectral weight in the 
first Brillouin zone as measured by ARPES for various 
cuprates. The spectral weight was integrated over an energy 
window of 30meV below the Fermi energy. Results are shown 
for a) Lai.6_ a; Ndo.4Sr a ;Cu04 (a; = 0.15) measured at 15 K 
(from Ref. M), b) I^-zSr^CuCU (x = 0.15) measured at 15 K 
(from Ref. m), and c) overdoped Pb-doped Bi2Sr2CaCu208+i 
(T c =70 K) measured at 20 K (from Ref. 0). 

In Fig. |^ we report results for an ordered array with 
a spacing between stripes of t = 4. In Fig. ^, we do the 
same for an array with a distribution of stripe spacings 
chosen randomly between 1 and 7, i.e., a flat distribution 
with mean I = 4. In both cases, we have fixed such 
that the average electron density per site is 1 — 5 where 
S = 1/21. {8 is known as the density of "doped holes.") 
This corresponds on the average to quarter filled stripes, 
such that for t = 4 we obtain S = 1/8. In these figures, 
we have chosen Aw = 0.2, but the qualitative character 
of the distribution is not highly sensitive to this choice. 
The data for the disordered array is for a given realization 
of the stripe distribution, but the system is large enough 
that the results are self averaging in the sense that the 
corresponding figures look similar to the eye for different 
realizations. 



3 



7T 



71 
T 



71 




10 20 30 JO 50 60 70 80 



X 



7T 

T 



0.S5 

. 925 
. 9 

0.875 
0.85 

. 825 
0.8 

0.775 




71 V 



TT 




20 



i 



60 



80 



FIG. 2. a) The low-energy spectral weight in the first 
quadrant of the Brillouin zone integrated over an energy in- 
terval Aw = 0.2 below [i for a disordered stripe array with 
mean spacing I = 4 and m = 1 at 5 = 1/8 doping. The sys- 
tem size is 320 x 320 and the amplitude color is coded from 
blue through green to red, as in Fig. |l| b) Typical electron 
density profile of the array as a function of x (perpendicular to 
the stripe axis). While some of the dips correspond to stripe 
positions, others, such as the one at x » 53, correspond to 
bi-stripes. Thus, in comparison to Fig. ^, below, the density 
of dips is smaller although the density of stripes is the same. 



Clearly, a number of salient features of the low-energy 
spectral weight of these various stripe arrays are remi- 
niscent of those measured in ARPES. In particular the 
disordered array, Fig. |^, looks strikingly like the corre- 
sponding spectrum in LSCO. In both cases, most of the 
spectral weight is in a broad, flat region close to the (tt, 0) 
point with a fainter image near (0, 7r). For the ordered 
array this exhausts the low-energy spectral weight, but 
for the disordered array, as in the experiments, a small 
amount of spectral weight lies along what would have 
been the band-structure Fermi surface in the nodal re- 
gion. (Of course, in order to compare the results more 
closely with the ARPES measurements we should sym- 
metrize these results around k y — k x to allow for stripes 
running in both directions in different domains.) 



FIG. 3. Same as in Fig. || but for an ordered stripe array 
with i = 4, S = 1/8 and Auj = 0.2. 



These results demonstrate that two-dimensional struc- 
tures, including an apparent Fermi surface, can appear in 
the ARPES data of a system which is dynamically one- 
dimensional. It is also noteworthy that y directed stripes 
not only lead to substantial low energy spectral weight 
in the anti- nodal region near (0, tt), as naively expected, 
but also in the reciprocal region, near (jt, 0), as well. 
(This observation is important in systems with macro- 
scopic stripe orientational order, including, potentially, 
YBCO [^2|.) We also note that for the ordered array, 
there are oscillations in the spectral weight as a function 
of k x . It is apparent from Fig. || that these oscillations re- 
flect the periodicity of the array (27r/4 = 7r/2 in this case. 
We have checked jl3| that for longer period arrays, for 
example with I = 6, the period of these oscillations shifts 
accordingly.) That such oscillations have been observed 
in LNSCO M further corroborates the stripe interpreta- 
tion of the results. 
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a large number of stripes, so we have purposely reduced 
the system size to L x = 80 and the energy window to 
\E\ < 2 for graphical clarity. We devote the next section 
to a closer examination of this band-structure. 



E 




FIG. 4. The squared magnitude of the Fourier transform of 
the field m(x, y) for a disordered array of size L x = 100, 000. 
This quantity is proportional to the spin-spin correlation func- 
tion measured by neutron scattering. 

We have also calculated the electronic density distri- 
butions of the disordered and ordered arrays and they 
are shown next to the spectral maps in Figs. [2] and |^. 
We see that in the ordered case, the density modulations 
have the same periodicity as the stripe potential, but in 
the disordered array, the mean spacing between troughs 
(or peaks) of the density is larger than the mean spacing 
between stripes. This reflects the fact that only a single 
density depression occurs where two stripes are close to- 
gether. To demonstrate this fact and make connection 
with neutron diffraction experiments, we have computed 
structure factors from the Fourier transforms of the den- 
sities of electronic charge (p) and z-spin component of the 
spin (S z ). For the ordered array of period 4 the charge 
signal is strongly peaked at k x = with a small satellite 
at k x = 7r/2 and the spin signal is peaked at k x = 37r/4 
and k x = 57r/4. For the disordered array the charge and 
spin peaks are shifted towards k x — and k x = n re- 
spectively In Fig. |J we show the Fourier transform of 
the autocorrelation function of the random field m. (We 
show the results for m rather than for S z , because it can 
be computed for very large arrays, where the result is 
self-averaging. However, we have verified that the two 
quantities give similar results jig] .) The fact that for 
the disordered array, the spin peaks are not only broad- 
ened, but are shifted towards the (tt, tt) point seems to 
be generic behavior for disordered arrays. 

Since k y is still a good quantum number, we can still 
talk about a (one-dimensional) band-structure of the dis- 
ordered array, as shown in Fig. [5]. However, as there are 
precisely 2L X bands, this figure becomes rather dense for 



FIG. 5. Band-structure of a disordered array with mean 
stripe spacing I = 4 and m = 1. p is the T = chemical 
potential for a density of doped holes 8 = 1/8. The system 
size is 80x320 

Because for each value of k y , the effective Hamilto- 
nian is that of a one-dimensional disordered system in 
the x direction, we know a priori that all the states are 
exponentially localized transverse to the stripes. How- 
ever, for the states that lie within the Mott gap, \E\ < 
m = 1, one might expect the states to be highly local- 
ized in the neighborhood of one or two stripes, while for 
\E\ > m = 1, the states are more likely extended scat- 
tering states, that are only weakly localized. That this 
expectation is realized can be seen in at least three ways. 
It is clear from Fig. ||, that there is a dense set of states 
at E > m = 1, whereas for E < m = 1 there is a spider 
web of identifiable one-dimensional bands; the discrete- 
ness of these bands is a reflection of their large degree of 
localization in the transverse direction. By plotting the 
states in real space along x we have confirmed that all 
the states in the gap decay exponentially in a few lattice 
constants, whereas the scattering states are essentially 
delocalized. Finally, we have computed the participation 
ratio, P a ee J2r\^a{r)\ 4 ~ (2£)~\ where £ is the local- 
ization length. For a system of size 320 x 320, we find Ps 
for the states with \E\ > m of order P ~ 1/100 (which 
might get even smaller for larger system sizes), while for 
\E\ < to, typically P « 0.1 — 0.25. For moderate doping 
5 the Fermi energy always lies inside the set of highly 
localized bands. (See for example, Fig. || where the value 
of the Fermi energy corresponding to S — 1/8 is indi- 
cated.) It is clear that all of the spectral weight shown in 
Fig. H comes from these localized, and hence dynamically 
entirely one-dimensional states. 
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IV. ISOLATED STRIPES AND BI-STRIPES 

In order to get a better understanding of the nature 
of the localized states and how they give rise to the 
very characteristic two-dimensional distribution of spec- 
tral weight, we look at some simpler problems consisting 
of a single or a pair of proximate domain walls. 

Consider the situation in which the defects, i.e., one or 
more stripes, are confined to a region about x = 0. It is 
straightforward to explicitly write the wavefunctions for 
fixed k y in the asymptotic regions to the left and right 
of the defect; the eigenvalue problem is then solved by 
matching these solutions across the defect region. For 
\E\ < to, the asymptotic states are exponentially falling 
functions of x. For instance, to the right of the defect, 
the Hamiltonian (|l|) is diagonalized by 



0k v (x,y) = e 



— P - Kx e i (<i x + k yy) 



(5) 



where the quantities k, q, and a are the implicit functions 
of k y and the energy E which (for a = 1) satisfy 



= cos(g) cosIi(k) + cos(fe„) , 
E 2 = to 2 - 4sinh 2 (K) sin 2 (g) , 
e ia = [E + 2i sinh(K) sin(g)]/m 



(6) 
(7) 
(8) 



It follows from Eqs. §| that if (§) is a solution then so 
is (||) with q — > — q, a — > —a and the same energy. Thus 
the total wavefunction (x, y) is a linear superposition 
of the two. To the left of the defect the same equations 
hold after substituting k — ► — k and, in the case of an 
anti-phase defect, m — > —to. A fourth equation, which 
depends on the nature of the defect, is obtained by in- 
tegrating the Schroedinger equation across the defect re- 
gion and defines the eigenvalue problem, E — E(k y ). (If 
the defect region is sufficiently broad, there may be mul- 
tiple solutions to the eigenvalue equation, corresponding 
to multiple mid-gap bands.) 

Much of the interesting physics is already implicit in 
these relations. The distribution of spectral weight in 
the two dimensional BZ associated with a state of given 
crystal momentum k y is given by \^>k (k x )\ 2 , where \& is 
the Fourier transform of $. Clearly, this weight is peaked 
near k x = q(ky) in a region of width Ak x rs n(k y ). For 
small K, Eq. (ph implies that q(k y ) simply traces out the 
underlying "diamond" Fermi surface given by the unper- 
turbed band-structure. 

To be explicit, let us consider the case of a single site- 
centered anti-phase domain wall, (Hit ■ ITIT) of the 
type which appears in Eq. In this case, the eigen- 

value equation is readily derived 



E(k y ) = ±2tanh 2 («) cos(fcy). 



(9) 



The resulting band-structure and spectral weight distri- 
bution in the BZ is shown in Figs. || and ^. Note that the 



mid-gap band has energy zero precisely where k y = ir/2, 
so that for a hole-doped stripe, for which the Fermi en- 
ergy will certainly lie at negative energies, there will never 
be low-energy spectral weight at the nodal point. A 
generic feature of the spectral weight distribution that 
is apparent in the figure is the fact that the negative en- 
ergy states have more weight in the lower-half of the first 
quadrant of the BZ than in the upper-half. This feature 
is seen very clearly also in the disordered array problem 
(Fig- 1). 
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FIG. 6. Dispersion E(k y ) of the bound states of a single 
anti-phase domain wall with m = 1. 
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FIG. 7. Integrated spectral weight over all negative (a) and 
positive (b) energy bound states of a single anti-phase domain 
wall with m = 1. The apparent discontinuity across k y — n/2 
is due to the band crossing. 

It is particularly illuminating to examine this result in 
the small to limit. The four implicit equations can be 
solved to leading order in m to obtain (for < k y < n) 



q{k y ) 
n{k y ) 

E(ky) 



k 



y i 



2sin(A;j / ) ' 
±2K 2 (k y ) cos(k y 



(10) 
(11) 
(12) 



As advertised, q(k y ) precisely traces out the diamond 
Fermi surface. n~ 1 (k y ) ~ VF(k y )/m is the extent of the 
wavefunction in real space, implying that |4'fc !J (fc x )| 2 is 
quite sharply peaked near k x = q(k y ) in the nodal re- 
gions, where the Fermi velocity, vf, is large, and is much 
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more diffused in the anti-nodal regions, where vf — > 0. 
(Where vp = 0, the expression for k breaks down and 
we find k ~ v 7 ^-) These features are, to a large extent, 
independent of the nature of the defect, so long as m is 
not too large. It is also worth noting that the energy 
width of the mid-gap band is of order m 2 , although the 
gap itself is 2m, and hence much larger. 

The single anti-phase domain wall solution captures 
most of the general features of the array of domain walls: 
The spectral weight is concentrated about the diamond 
Fermi surface, it is largest in a broad anti-nodal region, 
and is larger in the lower-half than in the upper-half of 
the BZ quadrant. Indeed, the low-energy spectral distri- 
bution of the ordered stripe arrays is, to first approxima- 
tion, simply the sum of contributions of isolated stripes, 
even when the stripe spacing is only four lattice con- 
stants. 

However, the disordered stripe array is not simply a su- 
perposition of single wall states, especially in the nodal 
region. In fact, while it is clear from the band-structure 
of the disordered array shown in Fig. [|, that a large frac- 
tion of the mid-gap bands are very similar to the mid-gap 
band of a single isolated stripe (Fig. ||), there are other 
bands, especially those that dominate the spectrum at 
k y = tt/2 and E near fi, that look very different. 

Guided by the observation that the nodal weight is 
missing for an ordered array, it is reasonable to assume 
that this feature is related to local configurations with 
two or more domain walls in close proximity. Consider 
therefore the problem of two anti-phase domain walls in 
close proximity. Clearly, the statement of proximity is 
related to the decay length £ = 1/k of the single do- 
main wall problem. If the distance between the walls is 
greater than £, the walls can be considered independent, 
while for smaller distances they interfere. For m = 1 we 
find 1/k « 2 and only weakly dependent on k y , implying 
that domain walls separated by more then two sites are 
roughly independent. Fig. ^ shows the energy spectrum 
for the bound states of two anti-phase domain walls sepa- 
rated by one site. The spectrum is gapped and it is clear 
that doping such a system can give rise to low-energy 
spectral weight around k y — tt/2. In Fig. || we repeat the 
same calculation for two proximate stripes. 

Particularly interesting is the distribution of spectral 
weight in momentum space for this system as shown in 
Figure [h]. As with the single wall problem, the weight 
hovers around the band-structure Fermi surface, but with 
the distinct feature that the weight of the negative en- 
ergy states is shifted towards the zone center whereas 
the positive energy states are shifted towards the (n, n) 
point. We find that the magnitude of this shift increases 
with m. Intuitively it is of course not surprising that the 
negative energy states are shifted towards the zone center 
where the tight-binding energies are lower and conversely 
for the positive energy states. 
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FIG. 8. Spectrum E(k y ) of the bound states of an effective 
in-phase domain wall consisting of two one site anti-phase 
domain walls one site apart (Till ■ T ' ITTI) with m = 1. 
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FIG. 9. Spectrum E(k y ) of the bound states of a two site 
wide in-phase domain wall (TITI ' ' Till) with rn = 1. 
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FIG. 10. Integrated spectral weight over all negative (a) 
and positive (b) energy states of an in-phase domain wall 
consisting of two one site anti-phase domain walls one site 
apart with m = 1. 



V. INCLUDING INTERACTIONS 

So far we have been concerned with the distribution of 
spectral weight in the non-interacting model defined in 
Eq. (Q) when holes are introduced into its single-particle 
levels. Including the effects of interactions between elec- 
trons in a doped stripe is a natural extension of this 
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model. In the following we consider this issue for the 
case of the disordered stripe array. 

To the extent that the localized bands shown in Fig. || 
can be considered as independent we may calculate the 
low-energy and long-wavelength spectral response of the 
system as a sum over interacting Luttinger liquid contri- 
butions Aid{k,ui) 



A(k,uj) = \^kAkx)\ 2 A ld [k y - fc F ,^; 7 (fc F )] . (13) 



Here the sum runs over the localized bands that cross 
the Fermi energy. Every crossing point is a Fermi point 
kp = kp ■ y for the one-dimensional channel defined by 
the band. Each channel is also characterized by the k x 
profile of its wavefunction ^k F {k x ) and its Luttinger pa- 
rameters K c%s and velocities w CjS which are collectively 
denoted as "/(kp) in (O). 



The evaluation of the spectral function ( |13|) simplifies 
in the small m limit where the spectral weight is con- 
centrated (within a width of order m) along the diago- 
nal connecting the anti-nodal points and is roughly con- 
stant. The total spectral function is then proportional to 
A\d{k x + k y — 7r, uj; 7), where we have assumed that 7 is 
the same for all bands. As we see, in this case, the spec- 
tral function exhibits one-dimensional behavior over the 
entire BZ. This description is valid as long as the stripe 
array is sufficiently dilute such that the overlaps between 
different localized bands are small. 

For intermediate values of m performing the convolu- 
tion in jl^ ) is difficult owing to the complicated distri- 
bution of the low-energy spectral weight in momentum 
space. However, the spectral function continues to mimic 
one-dimensional behavior along directions that are not 
necessarily aligned with the stripes. To demonstrate this 
point we consider the spectral function along a cut in the 
BZ that runs in the nodal direction [from (0, 0) to (it, it)]. 
Since the zero-temperature A\d{k — kp, uj; 7) vanishes un- 
less |w| > min(w c , v s )\k — kp\ [|| it is clear that the impor- 
tant contribution to A(k, uj) along this cut, and at low en- 
ergies, comes from bands that cross the Fermi energy near 
k y = ir/2. Moreover, we have shown that these bands 
produce spectral weight that is concentrated along a line 
parallel to the anti-diagonal [(tt, 0) to (0,7r)]. Therefore, 
for our purpose, we may approximate the spectral dis- 
tribution by a Lorentzian distribution centered on this 
line with half-width at half-maximum of k. In Fig. |ll we 
present the spectral function which results from Eq. ( 1^) 
under the assumption of such weight distribution. 
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FIG. 11. a) The zero-temperate spectral function A(k,u>) 
calculated from Eq. ( |l3| ) assuming a Lorentzian distri- 
bution of spectral weight centered on a line parallel to 
the anti-diagonal. Here v c — v a = 1, K a = 1, 
7c = (K c + K' 1 - 2)/8 = 0.3 and k = 0.5 (which is what we 
find for a single stripe with m — 1.) k is measured along the 
BZ diagonal and relative to the Fermi point near (n/2,n /2), 
uj is measured relative to the chemical potential /j,. b) MDC 
at u; = 0, c) EDC at k = 0. 

The spectral function in Fig. |ll| shows a close resem- 
blance to one-dimensional spectral functions (compare 
for example with Ref . |^) . In particular it exhibits a strik- 
ing dichotomy between sharp MDCs (momentum distri- 
bution curves - cuts at constant uj) and broad EDCs (en- 
ergy distribution curves - cuts at constant k), which we 
believe is a telltale of electron fractionalization ||. It is 
interesting to note, however, that A(k, uj) in Fig. |ll| is 
smoother than the corresponding zero-temperature one- 
dimensional spectral function and does not diverge at 
uj = vk for k < like the latter. This divergency is 
smeared by the convolution in Eq. (13) and the result 
is a spectral function that resembles A\d(k,uj) but at a 
finite temperature. 

The present analysis is the first which combines 
the remarkable non-Fermi liquid features of the one- 
dimensional electron gas with the two-dimensional struc- 
ture of a Fermi surface. At a qualitative level, it justifies 
the analysis presented in Ref. in which the measured 
spectral functions of the high temperature superconduc- 
tors was analyzed in terms of the corresponding expres- 
sions for a one-dimensional Luttinger liquid. In particu- 
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lar, we have presented here a microscopic rationale, which 
was missing at the time of that earlier work, for the ex- 
istence of approximate one-dimensional dynamics in the 
nodal region of the Fermi surface. 

VI. FURTHER CONSEQUENCES 

In addition to the gross similarities between the exper- 
imentally observed spectral functions in LSCO, LNSCO 
and BSCCO and those obtained from simple stripe mod- 
els, there are a number of other aspects of the present 
results that are potentially relevant to the interpretation 
of experiment. 

A phenomenon that is fairly generally observed in 
ARPES measurements on high temperature supercon- 
ductors is that a sharp quasi-particle like peak appears 
in the spectrum in the anti-nodal regions as the temper- 
ature is reduced below T c . Recent experiments on un- 
twinned YBa 2 Cu307_,5(YBCO) crystals have found that 
the spectral-weight of this sharp feature exhibits strong 
anisotropy between the different anti-nodal regions of the 
BZ jl2]]. There are good empirical Jl4|,|l5| and theoretical 
Jl6| reasons to believe that the stripes in this material 
are partially oriented by the chain potential on the Cu-0 
planes. Thus, a large anisotropy in plane-related fea- 
tures is to be expected. Moreover, we have previously 
shown B that the emergence of such a peak is naturally 
explained as a dimensional crossover in a quasi-one di- 
mensional superconductor. However, the presence of low- 
energy spectral weight in the anti-nodal region transverse 
to the stripes seems surprising at first. Nevertheless, as 
noted above, stripes induce low-energy spectral weight in 
both anti-nodal regions, although with different spectral 
weights - this is a consequence of strong mixing between 
states separated by the antiferromagnetic wave-vector, 
(•7T, 7r). The observed anisotropy is reproduced naturally 
by the present model - see Fig. ||. 

One would generally expect electrons diffracting off an 
ordered stripe array to exhibit intensity modulations with 
the corresponding period, as indeed is apparent in the 
period n/2 structure seen in Fig. ||. We believe, in agree- 
ment with the interpretation proposed Zhou et al. [Q, 
that this effect has already been seen in ARPES spectra 
of LNSCO @. Since all previous observations of charge- 
density wave ordering are indirect (e.g. they detect the 
accompanying lattice distortion density wave), this ob- 
servation is potentially of more far reaching importance 
as a method for detecting charge order. 

Experiments on both electron and hole-doped cuprates 
indicate that the chemical potential remains within the 
Mott gap, and is a remarkably slowly varying function 
of dopant concentration; instead, new dopant induced 
states appear in the gap Jl7|-Jl9[. Of course, this behav- 
ior is characteristic of the sort of stripe models envisaged 
here, in which the density of stripes increases with dopant 



concentration so as to accommodate the added charges 
in mid-gap states. This is a form of "topological doping" 
pot analogous to soliton doping of polyacetylene pll . As 
is the case in polyacetylene, we expect a variety of ad- 
ditional aspects of topological doping to be manifest in 
experiments on the cuprates, including dopant induced 
mid-gap optical absorptio n p] ,p||, dopant induced infra- 
red active phonon modes J22[| , etc. 

Finally, the observation that the structure factor of a 
disordered stripe array exhibits incommensurate peaks 
that are not only broadened by disorder, but also shifted 
to smaller incommensurability (Fig. ^) may have other 
interesting implications. Inelastic neutron scattering ex- 
periments in LNSCO LSCO || and YBCO || 
reveal incommensurate magnetic peaks displaced from 
the antiferromagnetic wave vector, (-7T, 7r), by an amount 
which first grows roughly linearly with doping, but then 
tends to saturate (typically beyond 5 = 1/8) . This is gen- 
erally taken as evidence that the concentration stripes, 
likewise, first increases linearly and then saturates. Our 
present results suggest that, at least in part, the satura- 
tion of the incommensurability may reflect increasingly 
strong stripe fluctuations, rather than a sudden satura- 
tion of the stripe concentration. 
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